Draft version May 3, 2009 

Preprint typeset using 1^1^^ style cmulateapj v. 10/09/06 



HELIUM-CORE WHITE DWARFS IN THE GLOBULAR CLUSTER NGC 6397 

Rachel R. Strickler^'^, Adrienne M. Cool^, Jay Anderson^, Haldan N. Cohn^, Phyllis M. Lugger*, Aldo M. 

Serenelli'' 
May 3, 2009 

Draft version May 3, 2009 



o 
o 

(N 



< 

o 

6 

(N 
> 

cn 

^' 

o 

o\ 
o 



X 



ABSTRACT 

We present results of a study of the central regions of NGC 6397 using Hubble Space Telescope's 
Advanced Camera for Surveys (HST ACS), focusing on a group of 24 faint blue stars that form a 
sequence parallel to, but brighter than, the more populated sequence of carbon-oxygen white dwarfs 
(CO WDs). Using F625W, F435W, and F658N filters with the Wide Field Channel (WFC) we show 
that these stars, 18 of which are newly discovered, have magnitudes and colors consistent with those 
of helium-core white dwarfs (He WDs) with masses ^0.2 — O.3A^0. Their Ha — i?625 colors indicate 
that they have strong Ha absorption lines, which distinguishes them from cataclysmic variables in the 
cluster. The radial distribution of the He WDs is significantly more concentrated to the cluster center 
than that of either the CO WDs or the turnoff stars and most closely resembles that of the cluster's 
blue stragglers. Binary companions are required to explain the implied dynamical masses. We show 
that the companions cannot be main-sequence stars and are most likely heavy CO WDs. The number 
and photometric masses of the observed He WDs can be understood if ~1 — 5% of the main-sequence 
stars within the half-mass radius of the cluster have white dwarf companions with orbital periods in 
the range ^1 — 20 days at the time they reach the turnoff. In contrast to the CO WD sequence, the 
He WD sequence comes to an end at i?625 — 24.5, well above the magnitude limit of the observations. 
We explore the significance of this finding in the context of thick vs. thin hydrogen envelope models 
and compare our results to existing theoretical predictions. In addition, we find strong evidence that 
the vast majority of the CO WDs in NGC 6397 down to T^ff ~ 10,000 K are of the DA class. Finally, 
we use the CO WD sequence together with theoretical cooling models to measure a distance to the 
cluster of 2.34 ± 0.13 kpc, or (m - M)g25 = 12.33. 

Subject headings: binaries: general - globular clusters: individual (NGC 6397) - stars: imaging — 
stars: Population II - stellar dynamics ~ white dwarfs 



1. INTRODUCTION 

Hubble Space Telescope (HST) has opened a new 
era in the study of white dwarfs by making it pos- 
sible to study these stars within globular clusters for 
the first time. Early HST observations revealed hand- 
fuls of white dwarfs in several clusters; since then, 
extensive cooling sequences comprising hundreds of 
whit e dwarfs h ave been o bserved in nearby clusters— 
e.g..lRenzini et al. (199 61): ICool. Piotto. fc Kinil (I1996D: 
[Richer et al.' (1997); Zoccali et al.' (2001*): [Hansen et all 
12004, 2007); Calamida et al. (2008). The advantages of 
studying a population of stars all at the same, well-known 
distance can thus be brought to bear on the study of the 
final stage of evolution of low-mass stars. 

As the second-nearest globular cluster to the Su n, and 
one with a well-constrained reddening (Gratto n et al.l 
120031 ). NGC 6397 is particularly favorable for the study 
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of white dwarfs (WDs). Early HST Wide Field Plan- 
etary Camera 2 (WFPC2) observations revealed a nar- 
row sequence of stars whose position was consistent with 
carbon-oxygen WDs (CO WDs) of a single mas s in the 
range ^0 5-0.6.^^ (ICool. Piotto. fc Kind [l996l) . Most 
recently, [Hansen et al.l ( 20071 ) have measured the dis- 
tance and age of the clus ter by det e cting the bottom 
of the WD sequence, and Davis et al.l ()2008f) have found 
evidence for natal kicks among CO WDs in the cluster. 
Both studies make use of deep H ST/ACS imaging o f a 
field outside the half-mass radius (|Richer et al.ll2008l ). 

NGC 6397 is also the nearest apparently core-collapsed 
globular cluster (Trager, King fc Djorgovski 1995). It 
has a small, resolved c ore of uncert ain radius (1.5 — 
8" at 95% confidence— ISosinI (|1997D 1 surrounded by 
a power-law "cusp" region extending out to ~100" 
(jLauzeral et al.lll992t iLugger et anil995f ). A wide va- 
riety of objects whose origins are thought to be linked 
to its high central density and resulting high rate of 
stellar interactions have been discovered in this cen- 
tral c usp region. These include numerou s blue strag- 
glers ([Auriere , Lauzeral. fc Ortolam ' 1990 1), cataclysmic 
variables (C Vs) (Cool et al. 1995; Grindl av et 10120011 
ICohn et al. ' 2009), a quiescent low-mass X-ray binary 
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A recent addition to the growing list of stellar exotica 
in NGC 6397 is a set of faint blue stars thought to be 
low-mass white dwarfs with helium cores (He WDs) . The 
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first three of these were found serendipitqusly i n a search 
for CVs using HST/WFPC2 (jCool et al.lll998h . Dubbed 
"non-flickerers," they showed neither the short-term vari- 
abihty nor the UV excess typical of CVs. Three more 
candidates were found in a subsequent WFPC2 study 
that included Ha imaging; all six candidates were found 
to be Ha- faint, c onsistent with thei r proposed classifica- 
tion as He WDs (jTavlor et al.ll200lD . The one candidate 
for which a spectrum has been obtained (NF2) revealed 
a broad H/3 line th at confirmed its identity as a low-mass 
white dwarf (Ed monds et al.|[l999D . 

Helium-core WDs have the potential to provide new 
insight into a cluster's population of binary stars — stars 
that play a crit ical role in globular cluster dynamics 
(|Hut et al.lll99 2^. The six He WD candidates identi- 
fied in NGC 6397 as of 2001 are confined to the central 
~30" of the cluster, a distribution that stron gly suggests 
the presence of unseen binary companions (jCool et al.l 
Il99a iTavlor et al.|[200ll ). The spectrum of NF2, with 
its apparently Doppler-shifted H/3 line, also hints at the 
influence of a dark, compact companion (Edmon ds et al.l 
|1999| ). Double-degenerate He WD binaries are of addi- 
tional interest as potent ial progenitors of AM CVn stars 
(iNelemans et al.l 1200 IbD an d sources of gravitational ra- 
diation ( Benacauistall2006l ) . 

In the field, He WDs are commonly found in 
binary systems, typical ly in ultra-compact X-ray 
binari es (ULXBs; INelemans. Jonker. fc Steeghsl 
(|2006( l and references the re in), as companions to 
MSPs ()van Kerkwiik et all I2005D or in d ouble- 
WD systems (iMarsh. Dhillon. fc Duckl Il995l : 
iMaxted. Marsh, fc MoranI I2002T In globular clusters, 
eight He WDs have been identified other than those in 
NGC 6397. Three are in ultra-compact X -ray binaries — 
in NGC 6712 (lAnderson et al.l Il99l, NGC 6624 
(I Anderson et al.l Il993^ and M15 (jDieball et all 
l2005f). Three are c ompanions to MSPs — in 47 Tuc 
(lEdmonds et"^l200l. NG C 6752 (|Ferraro et al.ll200l 
and M4 (jSigurdsson et all I2003D . In M4 a second 
He WD has been foun d as the companion to a subdwarf 
B star (jo 'Toole et al.l l2006). while in 47 Tuc a second 
was i dentified recently in a spectroscopic study of blue 
stars (|Knigge et al.ll2008f ). 

From an evolutionary point of view, the connection to 
binarity follows from the need to prevent helium ignition 
at the tip of the red giant branch (RGB) — otherwise a 
CO WD will result, as for single stars. If a giant's en- 
velope is removed before helium ignition, then the de- 
generate core will be exposed and produce a low-mass 
white dwarf with a helium core. Removal of the enve- 
lope may be accomplished either via Roche-lob e overflow 
onto a pre-existing close binary companion (jWebbinkl 
Il975l ) or a collision between an RGB star and a com- 
pact object, leading to a common envelope phase fo l- 
lowed by envelope ejection (jDavies. Benz. fc Hills! 1 19911 ). 
Either way, the He WD is left in a binary. He WDs 
can also form from single stars if they can somehow 
lose a n unusually large amount of mass while on the 
RGB ()Castellani fc Castellanil Il993l) . However, as we 
will show, the radial distribution of the He WDs in NGC 
6397 strongly suggests that they have binary compan- 
ions, making such scenarios unlikely to apply in this case. 

Here we use the Advanced Camera for Surveys' Wide 
Field Channel (ACS/WFC) to study the central regions 



of NGC 6397 in more detail. Our primary focus in the 
present paper is on the He WDs. With the wider field of 
view, improved spatial resolution, and greater sensitivity 
of ACS/WFC relative to WFPC2, we probe a larger re- 
gion of the cluster to fainter magnitudes, and uncover a 
much larger population of He WDs than was previously 
known. Future papers will provide a complete analysis 
of the full data set and report on our search for opti- 
cal counterparts o f Chandra X-ray sources in these data 
(|Cohn et al.ll2009l in preparation). 

Our use of the Ha and i?625 filters also enables us to 
test for the presence or absence of an Ha absorption line 
in the atmospheres of the CO WDs in the WFC field and 
thus to classify them as either of the DA (hydrogen atmo- 
sphere) or non-DA type. The relative numbers of DA vs. 
non-DA WDs in globular clusters is of significant interest 
in view of recent findings that the incidence of non-DA 
stars app ears to be much lo wer in open clusters than in 
the field (jKalirai et al.ll2005h . and outstanding questions 
about the mechani sms that give rise to these differences 
(|Moehler fc Bonol ((2008) and references therein). The 
possibility that the deficit of non-DA stars in open clus- 
ters may be explai ned as a result of th e larger masses of 
open-cluster WDs ([Kalirai et al.ll2005l ) can be tested by 
observing globular-cluster WDs, whose masses are sim- 
ilar to those of field WDs. All WDs for which classifi- 
cations have been made in globular clusters have been 
found to be DAs (jMoehler et al.ll2000l . [200l . However 
the number classified to date is small (4-5 in each of two 
clusters) owing to the difficulty of spectroscopic observa- 
tions of such faint stars in crowded fields. 

Finally, we determine a distance to NGC 6397 using 
^300 of the best-measured CO WDs, together with a 
CO WD model cooling track kindly provided by and 
transforme d to the HST/ACS filt er system by P. Berg- 
eron (Berg eron. Wesemael. fc Bea uchamp 1995). 

We describe the data set, together with the photomet- 
ric analysis and calibration method, in §2. In §3 we 
identify He WD candidates using color-magnitude dia- 
grams (CMDs). We then consider other possible expla- 
nations for stars appearing in the parts of the CMDs oc- 
cupied by the He WD candidates. In §4 we compare the 
CO WDs to theoretical cooling tracks in a Ha — i?625 
vs. i?625 CMD to distinguish DA from non-DA WDs, 
and discuss the significance of our finding that most, 
and perhaps all, of the CO WDs are of the DA class. 
In §5 we present our distance determination. In §6, we 
analyze the spatial, mass, and age distributions of the 
He WDs, and use the results to glean information about 
the nature of their binary companions, the orbital peri- 
ods of their progenitor binaries, and the binary fraction 
in the cluster We also compare our results to predic- 
tions made bv lHansen. Kalogera. fc Rasiol ()2003l ) regard- 
ing the He WDs in NGC 6397. We summarize our results 
in §7. 

2. OBSERVATIONS AND PHOTOMETRIC ANALYSIS 

We obtained images using the ACS/WFC comprising 
a total of five short and five long exposures each in the 
F435W and F625W filters (hereafter B435 and i?625) and 
40 long exposures in the F658N filter (hereafter Ha). Ten 
single-orbit visits were made over the course of a year 
beginning 2004 July 16. During each visit, one short 
and one long exposure were taken in either -B435 (13 s 
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and 340 s) or Rq25 (10 s and 340 s), followed by four 
390 s or 395 s exposures in Ha. Due to the changing 
roll angle of the telescope over the course of the year, 
a mosaic constructed from all of the images in a given 
filter results in the pinwheel shape seen in Figure 1. The 
mosaic provides complete coverage out to a radius of ~l!5 
from the cluster center, partial coverage in two or more 
exposures out to ~2!5, and partial coverage in at least 
one exposure out to ^2'.9. F or compariso n, the half-mass 
radius of NGC 6397 is 2^33 (iHarrisH iggel. 

Photometry was carried out using the techniques de- 
veloped by J. Anderson that use an "effective" point- 
spread function (PSF) approach to model how the PSF 
varies with position on th e ACS/WFC detector. T hese 
methods are described byjAnderson & King (2000) and 
apphed to the WFC byjAnderson & King (2006). We 
constructed effective PSFs for each of the three filters 
from the present data set using a subset of bright, well- 
measured stars in each filter. 

To generate a star list that would be close to com- 
plete while also excluding most artifacts, we began by 
finding peaks in each individual FLT-format image. We 
then transformed all coordinates to a common reference 
frame and looked for matches within 0.5 pixels. In B435 
and i?625, objects that created significant peaks in at 
least three of the five images were considered potential 
stars. In Ha a larger number of matches was required, 
depending on the number of images in which the star 
could potentially have appeared (up to 40). Multiple 
passes of star- finding were made; during each pass, the 
stars found in the previous pass were removed, and addi- 
tional stars were searched for in the subtracted images. 
W hen finding fainter stars , we used the method described 
bv lAnderson et al\ ()2008l ) to avoid false detections due 
to the 'knotty' PSFs and diffraction spikes associated 
with bright stars. The final collated list contained 25,004 
stars. 

For stars with i?625 < 18.5, we measured magnitudes 
in each individual frame and then averaged together the 
results. This worked well for bright stars which make 
strong peaks in all images in which they appear. For 
stars bright enough to saturate in all the long expo- 
sures we adopted magnitudes derived from averages of 
the short exposures instead. So as to obtain a complete 
census of bright stars, magnitudes were measured even 
for stars that are saturated in the short exposures. Sat- 
uration in the short exposures reduces the accuracy of 
the magnitudes measured for stars with -R625 ^ 14 (see 
Fig. 2). 

For stars with -R625 > 18.5, we combined the informa- 
tion from all available frames to measure a single mag- 
nitude for a given star. This produced noticeably better 
results for faint stars that do not necessarily produce 
strong peaks in every frame. By incorporating informa- 
tion from frames in which a given star would otherwise 
not have been measured, this method eliminates the bias 
toward brighter magnitudes that would otherwise occur 
for faint stars. 

To calibrate the photometry we first selected one long- 
exposure FLT image in each filter and perform ed pixel 
area map co rrections. Next we used DAOPHOT (iStetsonI 
[I987l[l99l to measure magnitudes in O'.'5-radius aper- 
tures for 15-20 bright, unsaturated, relatively isolated 
stars in each filter. We converted the resulting aperture 



magnitudes to the VE GAMAG system usi ng the pre- 
scriptions provided by ISirianni et al.l (|2005l ) . We then 
found the offsets between these calibrated magnitudes 
and our instrumental magnitudes for the same stars. Ap- 
plying these offsets to all the instrumental magnitudes 
then yielded calibrated magnitudes for the entire data 
set on the VEGAMAG system. 

3. AN EXTENDED SEQUENCE OF HELIUM-CORE WHITE 

DWARFS 

In Figure 2, we show i?625 vs. i?435 - i?625 and i?625 
vs. Ha~R(i25 color-magnitude diagrams (CMDs) for the 
central regions of NGC 6397. Here we plot only stars that 
were detected in at least two frames in each filter; 15,842 
stars made this cut. We further required that, in each 
filter, a star's profile be reasonably well matched to the 
PSF at its locatio n using the "q'' qualit y-of-fi t param- 
e ter d escribed by lAnderson et al.l ()2008[ ) — see IStricklej 
(200_q) for details. This eliminated another 866 objects, 
leaving a final list of 14,976 stars that appear in these 
CMDs. These are the stars that we analyze in the re- 
mainder of this paper. For stars brighter than i?625 = 
18.5 we plot averages of magnitudes measured in indi- 
vidual frames. For all fainter stars we plot magnitudes 
derived from simultaneous fitting. 

The most prominent features in these CMDs are the 
main sequence (MS) and the RGB, both of which appear 
nearly vertical in Ha—Re25- NGC 6397's short blue hori- 
zontal branch, blue stragglers, and WD sequence can also 
be seen. In the right panel, these stars all appear to the 
right of the RGB or MS, consistent with the stronger 
Ha absorption lines associated with hotter stars. In 
— Rq25 vs. i?625, the majority of the stars imme- 
diately to the left of the MS are likely to be from the 
halo/bulge background, while stars to the right of the 
MS include both foregrou nd stars and MS-MS binaries 
(cf. iCool fc BoltolU 

Taking a closer look at the WDs in the left panel, we 
see that they form a densely populated sequence run- 
ning from i?625 ~ 22, ^435 - i?625 ~ down to i?625 
> 26, i?435 — i?625 ^ 1- These are the carbon-oxygen 
(CO) WDs associated with single-star evolution in the 
cluster — stars that have t urned off the main sequenc e 
within the pas t ~2 G yr (cf. ICool. Piotto. fc King! (|1996( l. 
[Hansen et al.l (|2007[ )) . In addition, one can see a group 
of objects that, for a given ,6435 — i?625 color, fie ~l-2 
magnitudes above the CO WD sequence. These stars 
are well separated in color from the halo stars, most of 
which have -6435 — i?625 > 1, and include the 6 stars 
previously identified as pro bable He WD s (see Fig. 3; 
ICool et al . (1998), Edmonds_eLaI| (|1999( ). iTavlor et all 
()2001h ). Many more stars can now be seen in this region 
of the diagram. These stars are the primary focus of the 
remainder of this paper. 

To distinguish potential He WDs from probable 
CO WDs we used the following procedure. We note 
at the outset that we do not expect to achieve a per- 
fect separation between the two populations. Our goal 
is to select a set of the best He WD candidates to study 
in detail while also retaining for consideration marginal 
candidates that are less readily distinguished from the 
CO WDs. To this end we began by interpolating a 
(I53A40 CO WD cooling track from the [Bergeron et aTl 
(|1995f ) models and fitting it to a subset of the WDs cho- 
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sen to be very likely CO WDs. Details of the fitting 
procedure are presented in §5. A best fit was obtained 
for a distance of 2.34 kpc; the resulting placement of the 
theoretical cooling curve is shown in the close-up of the 
WD region (Fig. 3). 

Next we selected all stars within ± 0.8 magnitudes in 
i?435 — -R625 color of the theoretical CO WD cooling track, 
a region which should enclose all possible WDs while 
excluding halo stars. For each star in this region we 
determined the offset between its -B435 — i?625 color and 
the -6435 — i?625 color of the 0.537^0 cooling track at the 
same i?625 magnitude. We then divided the stars into six 
bins: one bin for stars from i?625 = 18.9 — 23.5 and five 
bins 0.5 magnitudes wide from i?625 = 23.5 to 26.0. This 
division produced bins with ^70-100 stars each. Next 
we computed the standard deviation of the color offsets 
within each bin. To minimize the impact of outliers (e.g., 
He WDs) on our determination of a (which we wished to 
represent a for the CO WDs exclusively), we iterated 2- 
5 times (until convergence), eliminating stars with color 
offsets greater than 2.5cr at the end of each iteration. For 
a Gaussian distribution of CO WD colors, this choice 
would eliminate ~1 CO WD in each bin. Finally, we 
defined a preliminary list of the best He WD candidates 
to be stars that landed more than 3ct away from the 
cooling track on the red side. 

Two additional criteria were imposed before a star 
was put on the list of best He WD candidates. First, 
we required that it be brighter than i?625 = 24.5, since 
measurement uncertainties increase rapidly for magni- 
tudes fainter than this; the possibility that additional 
He WDs may lie at fainter magnitudes is considered in 
detail below. Second, we required that it be at least 
0.75 mag brighter than the 0.537^© cooling sequence 
to eliminate possible double CO- WDs; this last require- 
ment eliminated two additional stars. The 24 stars that 
met all these criteria are marked with large open cir- 

cles in Fig. 3; Th e six previously id entified candidates 

(|Cool et al.l 119981 iTavlor et all 1200 iD are among these 
and are shaded grey. These 24 He WD candidates are 
listed in order of increasing i?625 magnitude in Table 1. 
For each star we provide, in columns 1-8, an ID num- 
ber, the X and y coordinate in our i?625 mosaic, the R.A. 
and Dec. in J2000 coordinates, the radial offset from the 
cluster center in arcseconds, the i?625 magnitude, the 
i?435 — i?625 color, and the Ha — Re25 color. 

To identify probable CO WDs, we followed the same 
procedure but used a more stringent requirement in do- 
ing the sigma clipping {2.25a instead of 2.5cr), which 
yielded slightly smaller final values of a for each bin. We 
then defined CO WDs to be the 438 stars that lie within 
± 3(7 in any bin. These stars are marked as large black 
dots in Fig. 3. Stars that lie between the CO WDs and 
the best He WD candidates or that failed to meet the 
two additional criteria described above were retained as 
marginal He WD candidates. These 17 stars are marked 
with small open circles in Fig. 3 and listed in the lower 
portion of Table 1. 

Having defined regions in the CMD occupied primar- 
ily by CO WDs vs. He WD candidates, we wish to assess 
the extent to which poorly measured CO WDs might be 
included among either the best 24 He WD candidates or 
the 17 marginal ones. First we note that a close exam- 
ination of the i?435 — i?625 vs. i?625 CMD (Fig. 3, left 



panel) reveals 6 stars in the range i?625 = 22 — 24 that 
are offset to the left of the CO WD sequence by about 
the same amount that the marginal candidates are off- 
set to the right of the sequence. Thus it seems likely 
that the marginal group of 17 stars includes a few poorly 
measured CO WDs. By contrast, the best 24 candidates 
have no such counterparts to the left of the sequence and 
are likely to be a much cleaner set. 

To more thoroughly assess the level of contamination 
among the He WDs we used arti ficial star tests in a man- 
ner similar to that described bv lAnderson et aTl ([2008). 
We added artificial stars into the individual images, and 
then reanalyzed the images using procedures identical to 
those used to analyze the data originally (see §2). A total 
of ^^100,000 artificial stars were introduced with magni- 
tudes in the range i?625 = 20 — 28 and colors matching 
the O.53A^0 cooling sequence shown in the left panel 
of Fig. 3. We then used the artificial star results to con- 
struct 1000 simulated CO WD sequences by picking stars 
at random from the artificial stars that were recovered. 
For each simulated WD sequence, we extracted a number 
of stars equal to the number of observed CO WDs and 
with the same luminosity function, i.e., the same number 
of stars per 0.5 mag bin. In each realization, we counted 
how many of the artificial CO WDs had landed in the 
regions of the CMD previously defined to contain the 24 
best and 17 marginal He WD candidates, respectively. 

Among the 24 best candidates, we find that nearly all 
of the 18 stars with Rq25 < 23.5 are hkely to be genuine 
He WD candidates: the average number of CO WD inter- 
lopers in this region in the simulated CMDs was only 0.8 
± 0.8. Of the 6 stars with i?625 = 23.5 - 24.5, about half 
may be interlopers, judging from the simulated CMDs. 
Thus roughly 20 of the 24 best candidates cannot be 
explained as poorly measured stars; this accords with 
the results of a visual inspection which revealed no sign 
that these stars were unusually crowded in the images. 
Among the 17 marginal candidates, the simulated CMDs 
suggest that roughly half are likely to be poorly measured 
CO WDs. We counted 2.8 ± 1.7 interlopers among those 
with i?625 < 23.5 (vs. 9 observed) and 1.6 ± 1.2 inter- 
lopers in the range i?625 = 23.5 — 24.0 (vs. 4 observed). 
For i?625 > 24.5, however, the artificial CMDs produced 
interlopers in numbers comparable to the number of ob- 
served candidates. We conclude that few if any of the 
4 faintest marginal candidates are likely to be genuine 
He WD candidates. 

For the remainder of the paper we will focus primarily 
on the subset of the 24 best He WD candidates, the vast 
majority of which are well-measured stars, while keeping 
in mind that perhaps half of the brighter marginal can- 
didates could also be genuine He WD candidates. Before 
going on to analyze their properties, however, we con- 
sider other possible explanations for the presence of stars 
in the region of the CMD that they occupy. We note that 
there is little question that they are cluster members: a 
uniform spatial distribution, as would be expected for a 
group of foreground or background objects, is definitively 
ruled out by their strong degree of concentration toward 
the cluster center (see §6.1). Our preliminary look at 
the proper motions of these stars provides independent 
confirmation that nearly all are cluster members; a full 
analysis of the proper motions is, however, beyond the 
scope of the present paper. 
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We first consider the possibility that they may be cat- 
aclysmic variables (CVs), eleven of which are now known 
in NGC 6397 (see triangles in Fig. 3). Given that two 
of these CVs land in the region in the -6435 — i?625 
CMD occupied by the He WD candidates, it is nat- 
ural to ask whether some of the He WD candidates 
co uld instead be CVs; t his p ossibility was suggested 
by iTownslev fc BildstenI (|2002D as a possi ble explana- 
tion f or the three candidates identified by iTavlor et al.l 
(|2001f ). One argument against this possibility is that 
while all the known CVs were found as counterparts 
to X - ray sources in C handra imaging (jCrindlav et al.l 
[2OOII : ICohn et all 120091 ) . none of the 41 He WD candi- 
dates identified here have X-ray counterparts down to 
the limit reached with existing Chandr a imaging {hx 
< 3 X 1029 erg s'^) (|Cohn et al.l [2001 . Very few if 
any CVs in the field have luminositie s below this limit 
(jVerbunt et all (|1997[ ). lPretorius et all ( |2007[ ). and refer- 
ences therein). 

A second argument against the He WD candidates be- 
ing CVs comes from comparing the locations of known 
CVs with the He WD candidates in the Ha — i?625 vs. 
i?625 diagram (Fig. 3, right panel). Nine of the 11 CVs 
have Ha — i?625 colors significantly to the left of the 
main sequence. This implies excess emission through the 
Ha filter, an d in turn the lik e ly pre sence of an Ha emis- 
sion li ne (see iGrindlav et al.l (|1995| ) and lEdmonds et al.l 
(|1999( ) for the efficacy of this method in picking out 
emission-line stars). Even the two CVs that lie to the 
right of the MS in Ha — Re25 are still significantly Ha- 
bright compared to the WDs; comparison to the white 
dwarf sequence is appropriate for most of the fainter 
CVs, since their blue -8435 — i?625 colors imply that their 
broad-band spectra are dominated by the hot WD. Thus 
all eleven of the CVs show strong evidence of having 
Ha emission lines, which is a key sig nature of accreti on 
and nearly ubiquitous in CVs (e.g., iWilliamj (|l983i) ). 
By contrast, nearly all of the He WD candidates have 
Ha — i?625 colors consistent with strong Ha absorption 
lines, and very similar to those of the CO WDs in the 
cluster. Only one (ID=14) has an Ha — Re25 color close 
to any of the known CVs. 

A second possibility to consider is that the He WD 
candidates could instead be detached WD-MS binaries. 
Like CVs, such binaries would naturally occupy the re- 
gion between the MS and WD sequence ijLopez fc Cool I 
l2007f ). To explain the proximity of the candidates to 
the CO WD sequence in the ^435 — i?625 vs. i?625 CMD, 
however, would require that the MS companions all have 
extremely low masses. Combining t he observed CO WD 
sequence with model MS stars from iDotter et al1(|200l . 
we find that all but the brightest 4 He WD candidates 
could be explained as CO WD/ MS-star binaries o nly if 
the MS star had M < 0.15Mq. iCool et all (|1998D set a 
similar limit on the possible presence of MS companions 
for the three original bright candidates using multi-band 
WFPC2 data. That the Ha ~ i?625 colors of the He WD 
candidates closely resemble those of the CO WDs also 
implies that if any of the candidates were CO WDs then 
the MS companions must contribute a small fraction of 
the total light of the system. As there is no reason to 
think that binaries in NGC 6397 should favor such low- 
mass companions — and exchange interactions certainly 
would not — we conclude that it is unlikely that many of 



these stars are detached WD/MS binaries. 

4. HYDROGEN VS. HELIUM ATMOSPHERES 

Use of the Ha filter enables us to determine whether 
hydrogen is present or absent in the atmospheres of 
both the CO and He WDs we observe. In Fig. 3 we 
have plotted a DB (helium atmosphere) cooling sequence 
(thick solid line) along with the DA (hydrogen atmo- 
sphere) cooling sequence (thick dash line). We interpo- 
lated thes e tracks f o r an a ssumed mass of O.53A40 from 
the Bergeron et al.l (|1995D CO WD models. In the right 
panel, the same two tracks are also shown shifted 0.07 
magnitudes to the right (thin lines) for reasons that will 
be explained below. 

The DA tracks span the entire range of observed mag- 
nitudes, while the DB tracks beg in at T^ff = 30, OOOK 
(just below the so-called DB gap — iLiebert et al.l ([l986)). 
In -8435 — i?625 (left panel), the DA vs. DB tracks sep- 
arate by at most '^0.1 magnitude in color (near i?625 
~ 24) — too small a separation to definitively distinguish 
DA from DB WDs using the present measurements. In 
Ha — Re25 (right panel), the tracks are separated by >0.2 
magnitudes in the range i?625 = 22.5 — 25 and as much 
as ~0.4 magnitudes near i?625 = 24. Thus DA vs. DB 
WDs are more readily distinguished using this filter pair. 
Note that the DB tracks should be taken in this context 
to represent any non-DA-type WDs, i.e., any WD with 
no Ha line, rather than DBs exclusively. We also note 
that the location of the DB track to the right of the 
main sequence in the right panel of Fig. 3 is the result 
of these stars being significantly hotter than the stars on 
the lower main sequence. Because the Ha filter is cen- 
tered at A 6584 A ( i.e., at a longer wavelength than 
the i?625 filter which is centered at A ~ 6310 A), hotter 
stars will naturally produce lower ratios of Ha to i?625 
flux and correspondingly larger values of Ha — i?625 even 
in the absence of any lines. 

To determine which of the WDs in NGC 6397 are likely 
to be DA vs. non-DA WDs, we focused initially on the 
CO WDs (large solid dots in Fig. 3) and proceeded as 
follows. First we considered only stars brighter than -R625 
= 24.5; for these stars, the separation between the tracks 
in the Ha—Re25 vs. i?625 diagram is large enough relative 
to the spread in the measured colors of the CO WDs 
that we can distinguish DA from non-DA WDs (Fig. 3, 
right panel). An inspection of this diagram reveals that 
nearly all of the CO WDs follow the DA cooling track, 
while there is no sign of any concentration of CO WDs 
around the DB track. This strongly suggests that the 
vast majority of CO WDs in NGC 6397 are of the DA 
type. 

We note that while the shape of the observed CO WD 
sequence in the Ha — i?625 vs. i?625 CMD is very simi- 
lar to the shape of the theoretical CO WD track (thick 
dashed line), there is an offset between the color of the 
track and the mean color of the CO WDs. We estimate 
that the offset is ^0.07 magnitudes, and surmise that it 
may be the result of calibration uncertainties; shifted DB 
and DA tracks are shown as thin solid and dashed lines, 
respectively. The existence of this small offset does not 
affect our conclusions below regarding the relative abun- 
dance of DA vs. non-DA WDs in the cluster. Nor does 
it affect our measurement of the distance to the cluster 
(see §5) which relies only on the -B435 — i?625 vs. -R625 
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diagram. 

We can make a quantitative comparison between WDs 
in NGC 6397 and those in the field by focusing on a 
portion of the CMD f or which the DA to non-DA ratio 
in the field is known. iTremblav fc Bergeronl (|2008l ) find 
that in the range T^ff = 10,000 - 14,Q0Qif the ratio 
of DA to non-DA WDs is about 4:1. For the models 
we are using (and at the distance we find — see §5), this 
temperature range corresponds to a magnitude range of 
i?625 — 23.5 — 24.5. In this magnitude range we observe 
126 CO WDs, all but one of which lies closer to the DA 
track than to the DB track; this number increases to 
four if we compare instead to the shifted tracks. If the 
DA/non-DA ratio matched that in the field we should 
have seen ^--^25 non-DA WDs. The probability of finding 
4 when 25 are expected is 2.3x 10^*". Thus it is clear that 
the ratio of DA to non-DA WDs is considerably higher 
in NGC 6397 than in the field. Given the likelihood of a 
few poor measurements among the 126 CO WDs in the 
interval in question, it is possible that all the observed 
CO WDs in NGC 6397 have hydrogen atmospheres. 

This conclusion does not depend on our assumption 
that the CO WD mass is O.53A40, as the position of the 
tracks in the Ha — -R625 vs. i?625 diagram depends only 
very weakly on mass. For example, the O.51A^0 cooling 
track lies much less than 0.01 magnitudes to the right 
of the O.53A^0 track in Ha — i?625- Because the tracks 
are nearly vertical in Ha — i?625 vs. i?625, our conclu- 
sion is also largely independent of the adopted distance. 
Changing the distance modulus by up to 0.2 magnitudes 
would change the Ha — Re25 colors of the tracks by less 
than 0.01 magnitudes. 

These results demonstrate that the incidence of non- 
DA WDs is much lower in at least one globular clus- 
ter than in the field. That more than 96% of the 126 
CO WDs we observe in NGC 6397 appear to have hy- 
dr ogen atmosphe r es con fi rms and extends th e findings 
of iMoehler et all ()2004[ ). iKalirai et all (|2005i) found a 
similar preponderance of DA WDs in open clusters, and 
suggested that it could be the result of the higher masses 
of CO WDs in open clusters vs. the field. That the same 
phenomenon has now been seen to occur in a globular 
cluster, whose CO WDs have masses quite comparable 
to those in the field, makes it unlikely that mass could be 
the primary factor in determining the relative numbers 
of DA vs. non-DA stars. 

Our findings also illustrate the power of Ha imaging 
for efficiently classifying large numbers of WDs. This is 
particularly valuable in globular clusters in which WDs 
are faint and crowded by much brighter stars, making 
spectroscopic observations very challenging. In a single 
study we have increased the number of WD classifications 
in globular clusters by more than an order of magnitude. 
Ha imaging has been validated in its ability to pick out 
emission-line stars ( e.g., CVs) thr ough follow-up spectro- 
scopic observations (jGrindlav et al. 1995). Further test- 
ing of the method, in the context of absorption-line stars, 
would be valuable to gain greater confidence that it reli- 
ably distinguishes DA from non-DA WDs. 

Finally, we note that the He WDs also follow the DA 
cooling track (see Fig. 3, right panel), showing that they 
too have hydrogen atmospheres. Independent of any 
comparison to theoretical tracks, their distribution is 
well- matched to those of the CO WDs in Ha — i?625, and 



like them shows no tendency to follow the DB track. The 
only apparent difference is that the spread in Ha — i?625 
colors for the He WD candidates seems to be somewhat 
larger than that for the CO WDs. This could be due in 
part to their larger range in temperatures and masses at 
a given i?625 magnitude relative to CO WDs. We have 
explored the effect of these parameters, however, and find 
that they are unlikely to be able to fully explain the ap- 
parently larger spread; it is unclear at present what the 
cause might be, if in fact it is significant. 

5. THE DISTANCE TO NGC 6397 

The distance to NGC 6397 has been subject to con- 
siderable uncertainty, with reported values ranging fr om 
2.2 kpc (Anthonv- Twarog. Twarog fc Suntzeff l9921 to 
2.8 kpc (Reid 1998). Most recently. [Hansen et al. (.200^) 
have detected the bottom of the CO WD sequence in an 
off-center field in NGC 6397 and used it to determine a 
distance of 2.54 ± 0.07 kpc (converting from the quoted 
distance ir iodulus) and an age of 11.47±0.47 Gyr. Here 
we use the'Ber geron et al"] (fl995.) theoretical WD cooling 
models to determine a distance by matching the models 
to the observed WD sequence in Fig. 3. We assume a 
W p mass in the range 0.53 ± 0.02 Mq, as suggested 
bv lRenzini et al.l ([1996') based on the luminosities of the 
termination of the asymptotic giant branch (AGB), the 
RGB tip, the horizontal branch, and post-AGB stars. 
This mass range is also in accord with spectro scopic mass 
meas urements of CO WDs in NGC 6397 tM oehler et al.l 
HOOI, and encompasses the mass of O.51A^0 determined 
bv lHansen et al. (2007) (see their Fig. 19). In view of the 
good match between the DA cooling sequence and the lo- 
cation of the CO WDs in the Ha - i?625 CMD (see §4), 
we also assume that all the CO WDs are of the DA type. 

Starting with a 0.537^© DA WD model in the S435 - 
i?625 diagram, we tested distances from 2.20 — 2.80 kpc 
at 0.01 kpc intervals and found the best fit value as fol- 
lows. First we determined extinction values in both fil- 
ters assuming E{B — V) — 0.18. For stars with col- 
ors of CO WDs in NGC 6397, and adopting the redden- 
ing law of , Seaton, (1979), we obtained A435 = 0.75 and 
^625 = 0.48; these values are very weakly dependent on 
spectral type. We then shifted the theoretical coohng se- 
quence to account for the distance and for the extinction 
in both filters. Next we selected all stars within ±0.8 
magnitudes in i?435 — i?625 color of the theoretical se- 
quence and in the range i?625 = 18.9 — 26.0, excluding 
known CVs. We then computed the offset in ^435 — i?625 
color between each of these stars and the theoretical se- 
quence at the same i?625 magnitude. To select stars likely 
to be CO WDs, we computed the standard deviation of 
these color offsets (using the bins described in §3), iter- 
ating several times and removing stars with > 2a devi- 
ations at the end of each iteration. This fairly stringent 
requirement was imposed in order to reliably remove non- 
CO WDs from consideration. The best fit distance was 
taken to be the distance that yielded a mean i?435 — i?625 
color offset closest to zero. The best fit was achieved for 
a distance of 2.34 kpc, or an apparent distance modulus 
of {m — M)q25 — 12.33; the result is shown in both panels 
of Fig. 3 (long-dash lines). 

The main contributions to the uncertainty in the dis- 
tance come from uncertainties in the CO WD mass, the 
reddening, and the photometric calibration. Examining 
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the WD models, we find that a change of ± O.O2AI0 from 
a CO WD mass of O.53A40 corresponds to a change of 
~0.05 in the i?625 magnitude for the range of WD tem- 
peratures in question. This translates to a contribution 
to the uncertainty in the distance of ~0.05 k pc. The un - 
certainty in E{B — V) is 0.02 magnitudes f,Harrislll996f ). 
which translates to a uncertainties of 0.053 in ^525 and 
0.03 in — A(i25- Rerunning the model-fitting algo- 
rithm described above using E{B ~ V) = 0.18 ± 0.02, 
we find that the uncertainty in the reddening results in 
an uncertainty in the distance of '^0.05 kpc. Finally, a 
maximum uncertainty in the photometric calibration of 
up to 0.1 mag would correspond to a distance uncertainty 
of -0.11 kpc. 

Combining these contributions to the uncertainty in 
quadrature yields a final estimate for the uncertainty in 
the distance of 0.13 kpc. Thus, with an assumed CO WD 
mass of O.53AI0, we find a distance of 2.34 ± 0.13 kpc. 
We note that our final result is i n reasonable a g reeme nt 
with the distance determined by iHansen" "eiTaLl (jlOOT) if 
we adopt their value of O.^IMq for the mass of all but 
the very faintest CO WDs (see their Fig. 19). In this 
case our distance increases to 2.39 ± 0.13 kpc, which is 
to be compared with their value of 2.54 ± 0.07 kpc. 

6. DISCUSSION 

The sequence of 24 helium-core WD candidates identi- 
fied here quadruples the number known in NGC 6397. 
This is the first such extensive sequence of individ- 
ually identified He WDs found in a globular cluster. 
Their presence illustrates the complexity of WD se- 
que nces in GCs, wh i ch ha s been extensively explored 
by iHurlev fc Shara I ()2003l ). A large population of 
He WDs has also been proposed to explain the total 
numbers and range of colors of WDs observed in uj Cen 
ICalamida et al. 2008). Our ability to identify individ- 
ual He WD candidates in NGC 6397 presents a unique 
opportunity to study a large group of these rare stellar 
remnants all at the same, well-known distance. Given 
their connection to binary stars and/or stellar interac- 
tions, they provide a fresh window into the inner work- 
ings of a globular cluster. 

6.1. Photometric vs. dynamical masses of the He WDs 

We begin by examining the location of the He WD 
candidates in the -B435 — i?625 vs. i?625 CMD, compar- 
ing their magnitudes and colors to theoretical He WD 
models. These He WD models have been computed 
wit h the same code an d input physics as described 
by ISerenelli et al.' (200 ^, and their synthetic spectra 
according to TRohr mann et al.l (|2002l ). Fig. 4 shows 
model sequences for masses ranging from 0. 175-0. 45A10 
and Z=0.00 02 progenito rs, a metallicity appropriate for 
NGC 6397 (|Harridll996l ). Because magnitudes (and to 
a much lesser extent colors) depend somewhat on the 
thickness of the remnant's hydrogen layer (and cooling 
age much more so — see below), we have computed two 
sets of models characterized by different H-layer masses. 
Models with "thick" envelopes (left panel) have H-layer 
masses ranging from 5 x 10~^M-q for the least massive 
model down to 4 x 10~'^Mq for the 0A5Mq WD model. 
These models were obtained by assuming that binary 
evolution leads to st able mass transfer; co mputational 
details can be found in lSerenelli et al.l (|2002f ). In this case 



mass loss from the WD progenitor ends when the enve- 
lope mass becomes small enough that an extended struc- 
ture cannot be supported any longer and then shrinks 
within the Roche-lobe. Models thus computed represent 
an upper limit to the possible H-layer thickness that a 
He- WD model of given mass can have. In these mod- 
els, nuclear burning (mostly through the p-p chain) is 
a main source of energy along the WD cooling phase. 
On the other hand, "thin" envelope models (Fig. 4 right 
panel) were obtained by modifying the envelope structure 
of pre- WD models in such a way as to render hydrogen 
nuclear burning at the base of the envelopes negligible, 
while keeping the total WD mass fixed. In our models, 
this translates into "thin" H-layers with masses between 
1 X 10-^7^0 and 5 x 10'^ Mq. Our sets of "thick" and 
"thin" models represent the two extremes in terms of 
cooling timescales for He- WD stars. They also bracket 
reality in terms of mass-radius relations for these stars. 
Even thinner H-layers than those adopted in our models 
would not translate into noticeable smaller radii (fainter 
magnitudes). 

Comparison to either set of models shows that most 
of the 24 He WD candidates (large circles) have mag- 
nitudes and colors consistent with masses in the range 
A4 '^O.2-O.3Al0. This is in good agreement with the 
mass of ~O.25A40 determined for the one candidate for 
which a spectrum has been obtained (ID=2 in Table 1; 
lEdmonds et al.l (|1999D ). Comparing the two panels, we 
see that with thinner envelopes (right panel) the inferred 
masses shift to slightly lower values on average. This ef- 
fect is particu l arly p ronounced at the bright end (see also 
IHansen et al.l (|2003D ). such that the brightest candidate 
lands outside the region spanned by the thin-envelope 
models in the right panel. Within the measurement un- 
certainties, the remaining candidates are all within the 
region spanned by either set of models. This includes 
candid ates 4, 5, and 6, first identified by iTavlor et all 
(|2001h . whose ^555 — J814 colors measured from WFPC2 
images appeared noti ceably redder than O.2A^0 models 
(jSerenelli et al.l[200^ . We speculate that the discrep- 
ancy may be the result of a bias toward redder measured 
colors for stars near the detection limit in the WFPC2 
observations. In both panels of Fig. 3 the marginal can- 
didates (small open circles) lie along tracks predicted for 
O.35-O.45A^0 He WDs. 

An independent measure of the average mass of the 
He WD candidates as a group can be obtained by ex- 
amining their radial distribution within the cluster. The 
dynamical relaxation time at the half-mass radius (r/j = 
2!33) i s ~3 X 10^ y ears; at the center it is less than 10^ 
years (|Harrisiri996f ). As all the He WD candidates are 
well inside the half-mass radius (and half of the best 24 
are less than 30" from the center), their radial distribu- 
tion should refiect their masses. Indeed, a high degree of 
mass segregation in the cent ral regions of NGC 63 97 has 
been well documented (Kin g. Sosin. &: Coo]||1995l ). 

In Fig. 5 we plot the cumulative radial distribution of 
the 24 best He WD candidates (solid black line) along 
with that of several other populations in the cluster: 
CO WDs (short-dash line), turnoff stars (long-dash), 
blue stragglers (dash-dot), and the 17 marginal He WD 
candidates (solid grey line). We also plot the cumula- 
tive distribution that would be expected for a spatially 
miiform distribution of objects (dotted line). Applying 
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the Kolmogorov-Smirnov test to compare distributions, 
we find first of all that the 24 best He WDs have a 
2x 10~^% chance of being drawn from a spatially uniform 
distribution, which definitively rules out a population of 
foreground or background interlopers. They are also sig- 
nificantly more centrally concentrated than the CO WDs 
(0.02% probability of being drawn from the same popula- 
tion) and main-sequence turnoff stars (0.2% probability). 
This strongly suggests that, on average, their masses are 
in excess of the turnoff mass of ^0.8A4q. These compar- 
isons do not take account of incompleteness, except that 
we excluded CO WDs fainter than _R625 — 24.5 (the mag- 
nitude of faintest of the good He WD candidates) . The 
greater difficulty of detecting He WDs as compared to 
turnoff stars in the central regions would only make the 
difference in radial distribution more pronounced. Of all 
the groups to which we compare the He WDs, they most 
closely resemble the blue stragglers (44.1% probability), 
which have a range of mas ses '^O.8-2A^0 (jSaffer et al I 
[200llDeMarco etal .112001 . 

In contrast to the 24 best He WD candidates, the 17 
marginal candidates are not nearly as concentrated to the 
center as the blue stragglers (K-S probability = 0.16% of 
being drawn from the same population; see Fig. 5). They 
most closely resemble the CO WDs (94.3%), though are 
also compatible with turnoff stars (48.8%). The proba- 
bility that these 17 stars are drawn from the same pop- 
ulation as the 24 best candidates is only 1.0%. Given 
the significant degree of contamination that we expect 
among these marginal candidates (roughly half are likely 
to be poorly- measured CO WDs — see §3), the similar- 
ity of their radial distribution to that of the CO WDs 
is perhaps not surprising. Still, if the other half really 
were He WDs, and these had masses comparable to those 
of the 24 best candidates, one might expect a greater 
degree of central concentration. We conclude that ei- 
ther the degree of contamination is greater than we es- 
timated or that the He WDs that do exist among these 
marginal candidates may not have the large dynamical 
masses characteristic of the 24 best candidates. In view 
of the uncertain status of the marginal candidates, we 
hereafter focus our attention primarily on the 24 best 
candidates. 

6.2. Nature of the binary companions 

The mismatch between photometric masses deter- 
mined from the CMD and dynamical masses determined 
from the radial distribution for the 24 best He WD 
candidates can be understood if they have binary com- 
panions of sufficient mass. Such companions are ex- 
pected on evolutionary grounds (see § 1) and earlier ob- 
servations found evidence for them (ICool et al.l 119981 : 
lEdmonds et"all Il999l : [Tavlor et all l2001h . Here ' we use 
the measured dynamical vs. photometric masses of the 
He WDs to place new constraints on the nature of these 
companions. 

Given that the bulk of the He WD candidates have 
photometric masses in the range 0.2-0.3A^o, whereas the 
dynamical masses resemble those of the much more mas- 
sive blue stragglers (0.8 — 2A4q), the companions must 
make up this large mass deficit. We note that the blue 
straggler mass range found by De Marco et al. (2005) is 
based on six of the brightest blue stragglers in NGC 6397 
and thus may not be representative of the entire popu- 



lation of 24 blue stragglers in our sample. To estimate 
the typical blue straggler mass, we have performed cored 
power-law fits to th e radial distributi ons of turnoff stars 
and blue stragglers (jCohn et al.ll2009l ) . Taking these two 
groups to be in approximate thermal equilibrium results 
in a blue straggler mass of 1.2 ± O.IMq. Adopting this 
as the typical mass for the He WD binary systems sug- 
gests a companion mass of ~1A^0. This rules out main- 
sequence companions, as MS stars near the turnoff mass 
would dominate the light and cause the He WD bina- 
ries to be significantly brighter and redder than is ob- 
served. CO WDs with masses ^0.53MQ{\ike those cur- 
rently being produced in the cluster) are also ruled out, 
since they fall well short of the required companion mass. 
We have also compared the observed radial distribution 
of the He WDs to the expected distribution of 1.6A^0 
stars (i.e., the dynamical mass of a He WD binary con- 
taining a neutron star) and find that 1AA4q companions 
are ruled out, even when we allow for incompleteness in 
the He WD distribution. Thus we conclude that the re- 
quired typical mass of the unseen companions of about 
1A^0 is most compatible with heavy WDs. 

The non-detection in X-rays of any of the 41 He WD 
candidates further supports the view that the compan- 
ions are white dwarfs rather than neutron stars. NSs 
should have been spun up to become MSPs during the 
mass-transfer stage. For comparison, the 19 MSPs in 
47 Tu c have -lO^o - lO^i erg s'^ ()Bogdanov et all 
[20061 ): if any of the He WDs in NGC 6397 had compa- 
rable MSP companions they should have been seen in 
existing Chandra irnaging that reaches L^^ < 3 x 10^^ 
erg s-i (|Cohn et al.ll2009l) . 

Our finding that the companions are most likely to be 
white dwarfs i s in go od agreement with the conclusions of 
iHa nsen et all (|2003[ ). who argue on theoretical grounds 
that the companions to the He WDs in NGC 6397 should 
be WDs and not neutron stars (NSs). A key part of their 
argument is that mass transfer from a ~O.87V(0 turnoff 
star onto a NS would be stable, whereas onto a WD it 
would be unstable. Stable mass transfer would yield rel- 
atively wide binaries which would be vulnerable to rapid 
breakup in the dense environs of the center of NGC 6397. 
Unstable mass transfer onto a WD would instead yield 
tight binaries (following a period of common-envelope 
evolution) which would be more resilient. We note that 
the latter scenario would apply to mass transfer onto 
white dwarfs with mass less than ~1A^0, but not for 
more massive WDs — a point we will return to below. 

6.3. Formation rate of He WDs and implications for 
binary population 

We can use cooling ages from the He WD models 
to consider what the distribution of He WDs in the 
CMD implies about the formation rate of these stars in 
NGC 6397. This in turn can provide insight into the 
viability of different models and/or formation scenar- 
ios. The analysis is complicated by the fact that cool- 
ing ages are highly dependent upon the assumed mass 
of the residual hydrogen layer, with thicker layers lead- 
ing to slower rates of cooling. Because thick hydrogen 
layers support hydrogen burning, WD evolution along 
the cooling sequence for thick-envelope models is com- 
pletely dominated by nuclear burning on the very long 
timescales dictated by the pp-cycle. By contrast, thin H 
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envelopes render nuclear burning negligible so that for 
thin-envelope models, WD evolution is driven solely at 
the expense of the star's internal heat. The differing 
rates of cooling can be seen in Fig. 4 by comparing the 
left and right panels, where "+" symbols are used to in- 
dicate cooling ages. In the left panel, marked ages are 
1 — 13 Gyr in 1-Gyr intervals. In the right panel, they 
indicate ages of 0.001, 0.0025, 0.005, 0.0075, 0.01, 0.025, 
0.05, 0.075, 0.1, 0.25, 0.5, 0.75, 1.0, 1.25, 1.5, 1.75, 2.0, 
2.25, 2.5, and 2.75 Gyr. Thus, for example, the brightest 
of the three He WD candidates identified bv lTavlor et al.l 
(|2001l ) at i?625 — 22.5 has a cooHng age of ~ 8 Gyr ac- 
cording to the thick-envelope models (left panel) vs. an 
age of only 0.1-0.25 Gyr according to the thin-envelope 
models (right panel) . We caution that the very youngest 
ages shown for the thin-envelope models (on the order of 
1 Myr) are uncertain by a factor of a few due to the fast 
initial rate of evolution and the somewhat ad-hoc nature 
of the models. 

Here we examine the implications of thin vs. thick en- 
velopes in turn. It should be emphasized that these mod- 
els represent only two possible sets of H envelope masses. 
However, they represent the two opposite extremes both 
in terms of cooling ages and mass-radius relations and 
thus should bracket reality. 

Considering the thin-envelope models first (Fig. 4, 
right panel), we see that the 24 best He WD candidates 
have implied cooling ages ranging from <1 Myr for the 
brightest candidate to about 1 Gyr for the faintest. The 
implied formation rate is then '--^24 Gyr~^. This is higher 
than the rate at which RGB stars and compact objects 
collide in the central regions of NGC 6397. Assuming 
then that the He WDs are instead the products of pre- 
existing binaries, the formation rate is determined by the 
fraction of stars that have a close white dwarf companion 
when they ascend the giant branch. The implied binary 
fraction can thus be determined by comparing the rate 
of formation of He WDs to the rate at which stars are 
turning off the main sequence in the same region. Count- 
ing stars above the turnoff at i?625 = 16.1, we find ~700 
subgiant and red giant stars. Considering the ^1.5 Gyr 
spent in these stages combined ( Pols et al.lll998[ ). this im- 
plies that stars are turning off the main sequence at a rate 
of about 470 Gyr^^ at the present epoch. Similar rates 
are implied by the number of horizontal-branch stars ob- 
served (48, which given a lifetime of ^^0.1 Gyr implies a 
rate of 480 Gyr^^) and by the number of CO WDs with 
ages of 1 Gyr or less (^480 after correcting for incom- 
pleteness). Thus the implied binary fraction is 24/470 
or about 5%. This is a small enough fraction that it is 
not surprising that the horizontal-branch and CO WD 
populations do not appear to be depleted by losses to 
He WD formation. 

The ~5% binary fraction required to account for the 
observed He WDs in the context of thin-envelope models 
represents a lower limit on the total binary population 
in the central regions of the cluster, as it accounts only 
for the particular class of binary that yields He WDs. 
The <5-7% binary fraction found amon g MS-MS bina- 
ries in the central region of the cluster bv lCool fc Boltoiil 
(j2002) would be in addition to the He WD progenitor 
binaries. The total will be less than the sum if some 
MS-MS binaries are themselves progenitors of the bina- 
ries that lead to the He WDs, i.e., if a CO WD replaces 



one of the MS stars through an exchange collision be- 
fore the other MS star turns off the main sequence (or 
at least before it climbs far enough up the giant branch 
to fill its Roche lobe). Such exchanges are not uncom- 
mon in N-bod y simulations of high-density clusters (e.g., 
Ilvanova et all (|2006f l). 

Interpreting the He WDs in the context of thin H- 
envelope models implies a binary fraction that is quite 
plausible. However, it does present one problem. The 
implied cooling ages of the brightest four candidates are 
very young, in the range of only ~l-2 Myr. It would 
be surprising to catch any He WDs so early in their 
evolution, even considering the factor of a few uncer- 
tainty in the coohng ages of such young objects. These 
stars cannot just be a small number of interlopers unre- 
lated to the cluster. Three have already been sho wn to 
be proper-motion members (jCool fc BoltonI [20021 ) . and 
our preliminary analysis shows that the fourth is as well. 
Such young ages would imply an implausibly high rate of 
formation in the very recent past, nominally on the order 
of 2 Myr""'^. Even taking the uncertainty in their ages 
into account, the present-day rate of formation would be 
more than an order of magnitude greater than the av- 
erage rate over the past 1 Gyr. The current formation 
rate would also have to be higher than the total rate at 
which stars are turning off the main sequence at present 
(~0.5 Myr~^). Yet only a fraction of these turnoff stars 
can be members of binaries that will produce He WDs. 
Thus it seems highly unlikely that these bright stars can 
be explained as He WDs with thin H envelopes. 

Examining the cooling ages in the context of thick H- 
envelope models instead (Fig. 4, left panel), the He WDs 
would have ages ranging from <1 Gyr to ~ 13 Gyr. In 
this case, the formation rate problem for the brightest 
candidates disappears, owing to the much slower rate 
of cooling associated with thick envelopes. The implied 
rate of formation in the recent past is ~5 Gyr~^, easily 
compatible with a binary formation scenario, as it would 
require that only ^1% of the stars turning off the main 
sequence produce He WDs in the present epoch. This 
rate of formation is also low enough that collisions be- 
tween RGB stars and CO WDs could make a significant 
contribution to their production. The rate is also not too 
different from the long-term formation rate of ~2 Gyr~^ 
implied by seeing 24 He WD candidates with ages up 
to ^^13 Gyr. Indeed, a moderate increase in production 
rate would not be surprising considering the inevitable 
dynamical evolution of the cluster on this timescale. 

In summary, either thick or thin envelopes can explain 
the observations with reasonable implications for the for- 
mation rate of He WDs (^^2 — 24 Gyr^^) and the binary 
fraction in the cluster (^^1 — 5%). The one exception is 
that the brightest candidates cannot be readily explained 
in the context of thin envelopes, as they would require 
an exceedingly high rate of formation in the very recent 
past. 

Additional insight into the population of He WD pro- 
genitor binaries can be gleaned by examining the distri- 
bution of He WD masses, as inferred from their positions 
in the CMD. The mass of the degenerate core that will 
be left behind when a giant loses its envelope to a binary 
companion depends on how far up the RGB the giant 
has climbed before overflowing its Roche lobe. This in 
turn depends on the orbital period of the binary: longer 
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period systems, being wider, allow the giant (and the 
mass of its He core) to grow larger before the envelope is 
lost. The He WD masses we deduce from the CMDs are 
mostly rather small. All but 2 of the 24 best candidates 
have implied masses in the range ^0.2 — 0.3Mq. This in 
turn implies that the pr ogenitor binaries have periods in 
the range ~1 — 20 days (|Rappaport et al.lll995[ ). Longer- 
period binaries containing CO WDs and MS turnoff stars 
or giants appear to be significantly more rare, given our 
finding (see §3) that at most ^6 of the marginal can- 
didates with i?625 = 22 — 24 are likely to be genuine 
He WDs. 

6.4. Termination of the He WD sequence 

Irrespective of the thickness of H envelopes, the He WD 
sequence comes to a relatively abrupt end near i?625 ~24. 
The dearth of He WD candidates at fainter magnitudes is 
readily visible in the Fig. 4 CMDs (recall from §3 that our 
artificial star tests show that the four faintest marginal 
candidates are likely to be poorly measured CO WDs). 
By contrast, the number of observed CO WDs is contin- 
uing to rise below i?625 ~24. The lack of fainter He WDs 
cannot be the result of incomplete recovery of faint stars. 
Artificial star tests show that in the bins from i?625 = 
24.5-25.0 and 25.0-25.5, the completeness is still -60% 
and 50%, respectively (see Fig. 6, top panel). Thus it is 
unlikely that any significant number of fainter He WDs is 
present. When we correct for incompleteness (see Fig. 6, 
panels 2-4), the CO WD luminosity function continues to 
rise to the limit of our observations, whereas the number 
of He WDs is clearly in decline. 

To further examine the significance of the absence of 
convincing He WD candidates below -R625 ~ 24.5, we 
measured the offsets in color from the 0.53 M© CO WD 
cooling track for each of the CO WD and He WD can- 
didates. Fig. 7 shows color offset distributions for stars 
with i?625 = 19 - 24.5 (left panel) and i?625 = 24.5 - 26 
(right panel). The distribution of the brighter stars 
is clearly bimodal, with the He WD candidates form- 
ing a peak centered ~0.4 magnitudes redward of the 
CO WD sequence. For the fainter stars, the color spread 
is larger due to growing measurement uncertainties. Still, 
there is no sign of any excess of stars on the right side 
(A(_B435 — _R625) ^0.3) where He WDs would be expected 
to lie. We estimate that at most —10 He WDs (2f7 upper 
limit) could hide in the wings of the CO WD distribution 
in the magnitude range i?625 = 24.5-26. 

The significance of the termination of the He WD se- 
quence is different in the context of thin vs. thick H- 
envelope models. If the hydrogen envelopes arc thin, the 
sequence ends at ~1 Gyr, well below the age of the clus- 
ter. S uch a termination was predicted by Hanse n et aD 
(|2003f ). who suggested that within ~1 Gyr, binaries con- 
taining He WDs would either merge via gravitational- 
wave radiation or be broken up in exchange interactions. 
The He WDs, as the lowest-mass stars in most such ex- 
changes, would typically be removed from the binary and 
quickly redistribute themselves according to their single- 
star mass. While we find that the sequence extends about 
one magnitude fainter than predicted by [Hansen et al.l 
([2003) ( who suggested that n one fainter than those iden- 
tified by iTaylor et al.l (pOOl) would be found in deeper 
searches) , our finding agrees with their prediction within 
the large uncertainties in cooling age associated with the 



unknown H-envelope thickness. In particular, while their 
"moderate" H envelope niodels put the cooling age of 
the faintest iTavlor et al.l (|2001[ ) candidate at ~1 Gyr, 
the "thin" H-envelope models we use here result in more 
rapid cooling, such that the He WDs are a magnitude 
fainter at 1 Gyr. 

While our observ ations confirm the prediction by 
[Hansen et al.l ()2003[ ) that the sequence should end at 
~1 Gyr, there is still a puzzle. This prediction was made 
specifically for the cluster core, assuming a central mass 
density of 1.5x10^ Mq pc~^. Yet few, if any, of the 
He WD candi dates are act ually in the cluster core, which 
is very small. SosinI ()1997f ) found a core radius of 1.5 — 8" 
at 95% confidence, while our preliminary analysis of the 
present data set suggests a radius in the range ~5 — 7". 
Outside the cluster core, in the so-called "cusp" region, 
the density falls off approximately as (r/rc)~^. Given 
that a binary's lifetime to exchange is inversely propor- 
tional to the density, a binary that would last only 1 Gyr 
in the core would survive more than a Hubble time at just 
4rc. Depending on the core radius, between half and all 
the He WD candidates are outside 4rc, where the lifetime 
to exchange would be even longer. Thus it is not clear 
that the termination of the sequence can be explained 
in this way. N-body simulations that can more directly 
assess the longevity of binaries, without the simplifying 
assumption that they reside at a fixed radius, would be 
valuable. 

In the context of thick envelopes, as illustrated in the 
left panel of Fig. 4, the faintest of the He WD candi- 
dates would have coohng ages of ^13 Gyr. Very old 
He WDs could have evolved from main-sequence stars 
with masses of ^ ,2Mp), whose MS lifetimes are <1 Gyr 
(jPols et al.lll998[ ). Such stars have core i nasses <O.SM.p) 
as they are leaving the main sequence (iN elemans et al.l 
2001a). Using the present thick-envelope models, 
the total inferred ages would be uncomfortably high 
considering rec ent age determinations for the clu ster: 
12 ± 0.8 Gyr (lAnthonv-Twarog fc Twarog I l2000l) and 
11.47 ± 0.47 Gvr ("Hans en et al.ll2007[ ). This alone should 
not rule out thick envelopes, however. Given the strong 
dependence of cooling rates on H envelope mass, even 
slightly lower envelope masses would likely bring the ages 
into agreement. Moreove r, the argument agai nst thick 
envelopes put forward by ' Hansen et al.l (|2003f ) rests in 
large part on their finding that binaries in the core will 
either merge or be broken up after 1 Gyr. They also 
suggest that common-envelope (CE) evolution will lead 
to thin H envelopes. However, in view of the location 
of many of the binaries outside the cluster core, and the 
lack of direct evidence that CE evolution leaves behind 
thin envelopes, we conclude that thick envelopes are not 
ruled out. 

Another possibility is that the He WDs could be a mix- 
ture of systems, some with thick and others with thin H 
envelopes. If common-envelo pe evolution does l ead to 
He WDs with thin envelopes ([Hansen et a l.' 2003). then 
all systems for which a ~O.8A^0 giant overflows its Roche 
lobe onto a CO WD companion with M < IMq should 
have thin envelopes and cool off quickly. For giants with 
CO WD companions of larger mass, however, the mass 
transfer would be stable. If thick envelopes were the re- 
sult, then this subset of systems would cool much more 
slowly. The brightest systems in NGC 6397 could poten- 
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tially be explained as such higher-mass binaries. 

This hypothesis is appealing in that it could help make 
sense of two other ways in which the four brightest 
He WDs seem to set themselves apart from the remain- 
ing 20 He WD candidates. First, their apparently greater 
concentration toward the center (median radial offset of 
12" vs. 33" for the fainter 20) could be explained as the 
result of mass segregation. Second, the 1.6-magnitude 
gap between the four brightest He WDs and the remain- 
ing systems could be the result of two very different rates 
of cooling. Judging from the models shown in Fig. 4, the 
youngest thick-envelope stars would have ages similar to 
the oldest thin-envelope stars (^1 Gyr). He WDs formed 
via stable mass transfer should have relatively long peri- 
ods following the mass transfer episode: ~2 — 40 days for 
0.2 - O.SA^o He WDs (Rap paport et all [1995?) . These 
binaries would be more vulnerable to exchange interac- 
tions than post-common-envelope binaries which would 
have much shorter periods. Nevertheless, at the range 
of radii where the four bright He WDs lie ('^6-24"), sys- 
tems with 2-day periods could survive ~0.6-10 Gyr, while 
systems with 40-day periods could survive '^0.08-2 Gyr 
before t he He WD would be exchanged out of the bi- 
nary (cf. [Hansen et all ((200l '). We surmise that the few 
we see concentrated to the cluster center could be the 
remnant of a larger population of such systems. 

Finally, we note that the possibility that some He WDs 
could be torn away from their binary companions via ex- 
change interactions raises the question of where these 
stars would be at present, and whether a second popu- 
lation of He WDs might be observable in the cluster. It 
may be that such stars would acquire high enough recoil 
velocities during exchange interactions that they would 
simply be ejected from the cluster. Escape from the clus- 
ter would follow, for example, if a 0.2A4q He WD in a 
10-day binary with a 1A4q CO WD acquired a veloc- 
ity comparable to its orbital velocity of ~90 km s~^. If, 
however. He WDs are retained in the cluster following 
exchange interactions, their half-mass radius should be 
larger than that of the cluster as a whole {rh ~140"), due 
to their very low masses. 

To estimate how many such isolated He WDs would 
still lie within the field of view spanned by the present ob- 
servations, we first note that any thin-envelope He WDs 
older than 3 Gyr would be fainter than the detection limit 
(see Fig. 4, right panel). Scaling from the 24 He WDs 
with ages of 1 Gyr or less, we estimate that another ~50 
would have cooling ages in the range 1 — 3 Gyr. Using the 
radial distribution of O.25A^0 stars in th e Fokk er- Planck 
models of NGC 6397 developed bv'PuIil (fl996 l). and as- 
suming that all 50 of these He WDs have been exchanged 
out of binaries, we estimate that ^--^8 isolated He WDs 
would currently lie within a radius of 1.5'. This is just 
barely compatible with our finding above that up to ~10 
He WDs could hide in the broad wings of the CO WD 
distribution at very faint magnitudes in the CMD. 

Using the same approach we can ask how many isolated 
He WDs would be pre sent in the off-center field studied 
bv [Richer et~al] (|2008D . Assuming that He WDs could 
be detected with ages up to ~10 Gyr in that field, we 
estimate that ~7 isolated He WDs would be present in 
their full ACS/WFC field. Taking a look at the CMD for 
that field (see their Fig. 3), there is no clear sequence of 
He WDs, but there are nevertheless several stars in the 



region such stars would occupy. Since it is possible that 
these stars are unrelated to the cluster, we also inspected 
their proper-motion-cleaned CMD (their Fig. 5). In this 
CMD, only 2 stars are present in the He WD region. 
However, this CMD excludes roughly 40% of the area of 
the original field (the area not covered by the Ist-epoch 
WFPC2 data), such that we would predict the presence 
of ~4 isolated He WDs, consistent with what is seen. We 
conclude that current observations do not rule out the 
existence of a second population of very faint, isolated 
He WDs in NGC 6397. 

7. SUMMARY AND CONCLUSIONS 

We have identified a set of 24 faint blue, Ha-faint stars 
in the globular cluster NGC 6397 that are good candi- 
dates to be white dwarfs with helium cores. This is the 
first extended sequence of He WDs found in a globular 
cluster. The masses of the He WD candidates that we 
infer from comparisons to the colors and magnitudes of 
theoretical cooling models are in the range ^0.2— 0.3A^q. 
Dynamical masses inferred from the radial distribution 
are significantly higher, on the order of 1.2A^0, implying 
the presence of binary companions. Constraints from the 
observed colors of the He WDs and their lack of X-ray 
counterparts suggest that the companions are most likely 
to be heavy white dwarfs. 

The rate of formation of He WDs that can be inferred 
from their distribution in the color-magnitude diagram 
is highly dependent on their rate of cooling, which in 
turn depends on the uncertain thickness of their hydro- 
gen layers. We compare the observed systems to two set 
of models, one with thick and the other with thin H en- 
velopes, and infer that ~1 — 5% of stars currently turning 
off the main sequence eventually become He WDs. This 
can be understood if ~1 — 5% of main-sequence stars 
have or obtain CO WD companions by the time they be- 
gin evolving up the RGB. Initial binary periods of ~1 — 10 
days are required to produce He WDs with the masses 
we observe. We note that this binary fraction is a lower 
limit to the overall binary fraction in the cluster, as it ac- 
counts only for the subset of binaries that yields He WDs 
with CO WD companions. 

The He WD sequence ends about 1.5 magnitudes above 
the limiting magnitude of our study. The interpretation 
of this termination is different in the context of thick 
vs. thin H envelopes. In the case of thick envelopes, 
the faintest observed systems have ages comparable to 
the age of the cluster. In the case of thin envelopes, 
the faintest have ages of '--^1 Gyr. Neither the thick nor 
thin envelopes considered here provide a completely sat- 
isfactory explanation for the current observations. It 
does seem clear that the brightest four He WDs have 
thick H envelopes, as thin envelopes would imply very 
young ages that in turn would require an implausibly 
high rate of He WD formation at the present epoch. If 
the rest have thin envelopes then the termination of the 
He WD sequenc e mus t be explained. The suggestion 
bv [Hansen et al.l (|2003f ) that mergers and/or exchanges 
destroy the binaries after ~1 Gyr appears viable, and 
could perhaps be tested by wider-area searches for faint 
He WDs in off-center fields. If all He WDs have thick 
H envelopes, then the envelope thickness must be some- 
what lower than what we have assumed here in order for 
the oldest observed He WDs to have cooling ages com- 
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patible with the age of the cluster. 

Nearly aU of the CO WDs brighter than i?625 ^ 24.5 
{Tgff > 10,000k) have Ha — i?625 colors that strongly 
suggest that they have the broad Ha absorption lines 
characteristic of DA-type (hydrogen-atmosphere) WDs. 
This suggests that whatever mechanism is operating to 
produce DBs (or other non-DA type WDs) in the field is 
ineffective in NGC 6397. White dwarf mass is unlikely to 
be the primary factor determining the relative abundance 
of DA vs. non-DA WDs since the masses of CO WDs 
in NGC 6397 are similar to those of CO WDs in the 
field. The 126 stars in our sample of tentatively classified 
CO WDs represents the largest sample of WDs for which 
atmospheric classifications have been obtained in a single 
cluster. While further validation of the method is needed 
to ensure that it reliably distinguishes DA from non-DA 
WDs, Ha imaging should be a powerful tool to classify 
white dwarfs in star clusters. 

Finally, we have compared the observed CO WD se- 
quence to theoretical cooling models to measure a dis- 
tance to NGC 6397. With an assumed CO WD mass of 
0.537^0 and reddening of E{B - V) = 0.18, we obtain 
a distance of 2.34 ± 0.13 kpc. 
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Fig. 1. — Mosaic image of NGC 6397 in the S435 filter, con- 
structed from images taJjen on five visits spread over one year. 




Fk;. 2. — Color-magnitude diagrams for 14,976 stars in the cen- 
tral regions of NGC 6397. The white dwarfs on the bottom left 
side of the -R625 vs. B435 — R(i25 diagram (left panel) are the focus 
of this paper. In Ha — i?625 vs. -R625 these stars appear to the 
right of the main sequence, indicating strong Ha absorption lines. 
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Fig. 3.— White dwarfs in NGC 6397. The 24 stars that we 
identify as the best He WD candidates are marked with large open 
circles. These include 6 previously known candidates which arc 
shaded grey. Marginal Ifc WD candidates (another 17 stars) are 
marked with small open circles. Grey triangles mark the 11 known 
cataclysmic variables in the cluster. Stars selected as probable 
CO WDs are shown as large solid dots. All other stars are shown 
as small dots. The thick dashed and solid lines in both panels 
represent 0.53A^o DA and DB WD cooling tracks, respectively, 
shown at the best-fit distance found in §5. In the Ha — R62B 
vs. R625 CMD, the same pair of tracks is also shown shifted 0.07 
magnitudes to the right (thin lines) to better align the DA track 
with the observed WDs (see §4). 
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Fig. 4. — Same as Fig. 3 but with He WD cooling tracks 
from low-metallicity p rogenitors (Z=0.0002) overlaid (dotted lines) 
llSerenelli et al.l[2002l) . The left and right panels show models with 
thick and thin hydrogen envelopes, respectively (see §6.1). In both 
panels, tracks are for He WDs with masses of O.175A^0, O.2OjV(0, 
0.25A^Q, 0.35Mq, and O.45A^0 (top to bottom). In the left panel 
(thick envelopes), the "+" signs indicate ages from 1 — 13 Gyr 
in 1 Gyr intervals. In the right panel (thin envelopes), the "+" 
signs indicate ages of 0.001, 0.0025, 0.005, 0.0075, 0.01, 0.025, 0.05, 
0.075, 0.1, 0.25, 0.5, 0.75, 1.0, 1.25, 1.5, 1.75, 2.0, 2.25, 2.5, and 
2.75 Gyr. 
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Fig. 5. — Cumulative radial distributions of various populations 
of stars in NGC 6397 and for a hypothetical uniformly distributed 
population. These distributions are used for the Kolmogorov- 
Smirnov tests described in the text (see §6.1). The dark soUd line 
shows the distribution of the 24 best He WD candidates; the solid 
grey line is for the 17 marginal candidates. These can be compared 
with the distribution expected if the stars were unrelated to the 
cluster (dotted line) and with the distributions of 168 CO WDs 
(short-dash line), 1416 turnoff stars (long-dash line), and 24 blue 
stragglers (dash-dot line). 
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Fig. 6. — Results of artificial star tests. The top panel shows 
the fraction of stars recovered in artificial star tests for stars on 
the WD sequence as a function of i?625 magnitude. We used these 
completeness statistics to infer the numbers of CO WDs (panel 2) 
and He WDs (best 24 only; panel 3; all: panel 4) in the cluster 
down to il625 = 26. Shaded histograms indicate numbers observed; 
white histograms include the completeness corrections. Prom the 
two bottom panels we conclude that, in contrast to the CO WD 
sequence, the He WD sequence ends before the magnitude limit is 
reached. 
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Fig. 7. — Color distribution of CO WDs and He WD candidates, 
measured relative to the color of the CO WD cooling tra<;k shown 
in Fig. 3. The left panel shows stars with 18.9 < i?625 < 24.5; 
the right panel shows stars with 24.5 < i?625 < 26. Diagonal-line 
shading denotes CO WDs. The best He WD candidates axe shown 
in black and the marginal candidates in grey. Stars not assigned 
to cither class are white. Among the brighter stars (left panel) the 
distribution is bimodal, with the He WDs forming a distinct sec- 
ond peak. At fainter magnitudes (right panel), the broadening of 
the CO WD sequence due to increased measurement uncertainties 
makes it more diflBcult to discern He WDs. However, there is no 
sign of any excess of stars on the right side as compared to the left. 
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Helium-core White Dwarf Candidates. 
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Offset in arcscconds from the cluster center. For the center we adopted l?"* 40" 42?049, -53° 40' 
28'.'72 from[Sosin (1997), which corresponds to ( x, y) = (2900, 25 06) in the mosaic. ^ ID=1, 2, and 3 
correspond to NFl, NF2, and NF3 identified b y ICool et al.l | |1998|) , respectively. ^ ID=8, 10, and 15 
correspond to PC4, PC5, and PC6 identified bv lTavlor et al.l 1 I2OOID . respectively. 
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